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DEFINITIONS 


NAME 

Bubble Parameter 

Load Factor 
Void Fractton 
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Slip Ratio 


SYMBOL DEFINITION 


C' Empirical 'drag coefficient,' based on 

d experimental data, used to calculate 

slip force between gas and liquid. 

F Ratio of load (electrode) voltage 

to Induced voltage, V^/u^aB. 

a Ratio of volume occupied by gas to 

total volume, or flow cross-sectional 
area occupied by gas to total cross- 
sectional area. 

X Ratio of gas mass flow rate to total 

mass flow rate, m^/Grig + m^). 

K Ratio of gas to liquid velocity. 
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SUMMARY 


The revised amblent-temoerature NaK-nltrogen facility Is described. The 
maximum liquid flow rate and generator Inlet pressure are 10.9 kg/s (200 gpm) 
and 1.48 MPa absolute (200 pslg), respectively, compared with the previous 
values of 6 kg/s (110 gpm) and 0.72 MPa absolute (90 pslg). Satisfactory loop 
operation has been obtained, and new experiments with the second dlverglng-channel 
generator were completed. The principal experimental results were a higher 
power density for the same generator operating conditions, and an aoparent 
tendency for the efficiency to Improve more with Increasing quality at higher 
velocities than at lower velocities. 

An evaluation of an annular generator geometry Is presented. The advantages 
and disadvantages of the geometry are described, the equations developed, and 
solutions obtained for three cases—constant velocity and no armature reaction, 
laminar flow with no armature reaction, and armature reaction with constant 
velocity. Numerical examples show that: (1) the attainable terminal voltages 
appear to be very low, (2) flow reversal and large viscous loss occur at or 
below the desired power densities, and (3) armature reaction effects are Important 
and compensation techniques appear Impractical. Thus, this annular geometry does 
not appear attractive for either generator or pump operation. 

The Initial steps In the program to produce and evaluate liquid-metal foams 
are described. The future directions of the experimental generator program. 
Including foams, are discussed. 




I. INTRODUCTION 


The two-phase liquid-metal MHD generator program combines experimental and 
analytical studies In an effort to better understand the Dhyslcal processes 
occurring In two-phase generators and to establish the limits of generator 
performance. For the oast six years this program has been sponsored by the 
Power Branch of the Office of Naval Research. 

The generator experiments are conducted with an ambient-temperature 
NaK-N 2 test facility consisting of a two-phase mixer, MHD generator, magnet, and 
required auxiliary loop equipment. The original facility, used In all tests 
through July 1975, has been described elsewhere.^ The revised facility, com¬ 
pleted this year. Is described In Section II.1 of this report. Earlier experi¬ 
ments on generators with varying geometries have been performed and the results 
presented In previous annual reports. 2 ’^ 

Analytical models of the generator have been developed to aid In the basic 
understanding of the flow characteristics in the generator and to provide a 
(one-dimensional) computer code which Is capable of approximately predicting 
generator performance from the Independent Input parameters such as flow rate, 
magnetic field strength, and geometry. This computer code and the semi-empirical 
technique used to characterize the bubble parameter (churn-turbulent drag coefficient) 
used In the determination of the Interfacial slip ratio have been reported In 
detail. 2 ’ 5 

Previous generator experiments have used three different channels. Fig. 1.1. 

The first channel was designed for a silo ratio of unity as the actual slip ratio 
was not known, but higher slip ratios were observed experimentally and the liquid 

p 

velocity was not constant, decreasing along the generator. Following this a 

2 3 

constant-area channel was tested, and the velocity Increased along the channel. ’ 
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The second diverging channel, the third channel teted, has a smaller angle than 
the first channel, Fig. 1.1, to yield a close-to-constant liquid velocity at the 
design point. As described In the previous annual report,^ Improved generator 
performance was obtained with this channel, and the liquid velocity was almost 
constant. 

Previous experiments have been limited to liquid velocities of about 
12 m/s or less, smaller than anticipated for practical generators. In addition, 
recent experimental research has demonstrated that the slip loss due to the 
velocity difference between the gas and the liquid can be decreased by increasing 
the liquid velocity.^ For these reasons the facility has been revised to allow 
higher liquid velocities or flow rates and higher pressures. In this report the 
revised facility Is described along with the initial data obtained with the second 
diverging channel In the revised facility, Including some high-velocity runs. 

The test program, described In Section II, consists of the facility revisions, 
the NaK-nltrogen tests, and the Initial foam experiments. The results of the 
generator tests, presented In detail In Section III, Include a comparison of 
runs using the same generator with the original and with the revised facilities. 

An evaluation of an annular generator geometry Is in Section IV. The conclusions 
and suggestions for future work are given In Section V. 







II. TEST PROGRAM 


The experimental program is summarized In this section. The revised facility, 
designed for generator tests at higher liquid velocities and pressures. Is described 
In Section II.1. The program for the NaK-nltrogen generator Is In Section II.2, 
and the results are described In Section III. The Initial exDerlments using 
surface-active agents to produce a foam are summarized In Section II.3. 

II.1 Revised Ambient-Temperature NaK-NItrogen Facility 

The maximum liquid flow rate has been Increased from 6 kg/s to 10.9 kg/s 
(110 gpm to 200 gpm), thus allowing generator tests at velocities anticipated 
to be typical of practical generators. As Increasing the flow rate Increases 
the pressure difference across the generator, the available generator Inlet 
pressure was also Increased from about 0.72 to 1.48 MPa absolute (90 to 200 pslg). 

The original loop was constructed primarily of two-inch pipe and had 
several two-inch globe valves that caused severe flow restrictions. In addition 
the pressure tanks did not have a high-enough pressure rating for some of the 
proposed tests. Thus, the loop was completely dismantled and rebuilt, as shown 
In Fig. II.1 for the NaK flow circuit. The major changes In the NaK circuit were: 

1. Two new canned-rotor NaK pumps were added In series with the existing 
pump. The three pumps produce a gross pressure difference of 1.57 MPa (228 psl) 
at 10.9 kg/s (200 gpm). For experiments where the total pump capacity Is not 
required, one or two pumps may be bypassed. 

2. The existing seoaratlon and supply tanks were reolaced by a single, 
larger tank with a rated pressure of 1.31 MPa absolute (175 pslg). It has 
Internal baffles enabling It to separate twice the NaK and nitrogen flow rate 
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of the old tanks. A vortex suppressor consisting of a three-inch Inner pipe 
surrounded by a five-inch outer pipe so that the NaK flows up the annulus 
between the two pipes and down the Inner pipe Into the line to the pumps was 
Installed In the tank. 

3. The calibration tank was replaced by a larger one of 210 kg NaK 
capacity and a rated pressure of 2.14 MPa absolute (295 pslg). 

4. The 150 psl NaK filter housing was replaced by a 300 psl housing. 

5. The air-cooled heat exchanger was replaced with two water-cooled heat 
exchangers connected In parallel. This eliminated the pressure loss of the 
previous heat exchanger, and provides sufficient cooling capacity for year-round 
operation. (In fact. In warm weather the loop can now be made cool enough to 
condense moisture on the NaK pipes. Before, operation had to be curtailed In 
warm weather.) 

6. Most of the two-inch piping was eliminated, and the NaK flow control 
Is now provided by a three-inch air-operated valve rather than a two-inch manual 
globe valve. All non-control valves are ball valves, except for a four-inch 
Y-pattern globe valve In the Inlet line to the first pump. 

7. Most new components are carbon steel as stainless steel Is not 
required for a room-temperature NaK facility. 

8. All kerosene cooling was eliminated, and the magnet and load resistor 
are cooled by water. This prevents the possibility of a combined NaK and kerosene 
fire. 

9. Some minor NaK carryover with the nitrogen from the separator tank 
was observed. A three-inch commercial secondary separator has been ordered 
and will be Installed In the six-inch exhaust line. 
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Changes completed In the nitrogen circuit Include: 

1. The nitrogen purifier (NaK bubbler) tank was hydrostatically tested 
to Increase Its pressure rating. 

2. An electrically-operated valve to shut off the nitrogen supply In 
case of an emergency and a new pressure relief valve were Installed. 

3. Some pipe sections and valves were replaced to reduce the number of 
bends and fittings, and the pipe length, and thus reduce the pressure drop. 

4. A proportional controller for the nitrogen exhaust valves was added 
to hold the pressure constant In the separator tank. 

Extra nitrogen heater capacity will be required In the future for the higher 
flow rates. 

The modified loop allows two modes of operation that were not possible 
before. These are: 

1. Operation at up to 10.9 kg/s (200 gpm) NaK flow with a corresponding 
Increase In the nitrogen flow to determine the effect of liquid velocity on 
slip and generator efficiency. Extra capacity Is available and operation at up 
to 18 kg/s NaK flow Is possible. 

2. Variable back pressure at the generator exit, up to 1.14 MPa absolute 
(150 pslg) at reduced flow rates, to determine the effect of back pressure on 
generator performance. In this mode only one or two pumps may be used to avoid 
excessive pressure. Operation at low back pressures, down to about 0.12 MPa 
absolute (17 psla), Is also possible, and this allows larger pressure ratios 
across the generator. (The old facility was limited to pressures above about 
0.24 MPa absolute due to the high pressure drop between the separator tank and 
the suction side of the pump.) 
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Several safety features are Incorporated Into the modifications. The 
piping Is designed for 2.0 MPa absolute (275 pslg) and the separation/storage 
tank for 1.31 MPa absolute (175 pslg). Pressure switches stop or prevent from 
starting one, two, or all three pumps If pre-set pressures In the separation 
tank or high-pressure piping are exceeded. The separation and nitrogen purifier 
tanks are protected by safety valves, as before, but with higher pressure ratings. 
Emergency switches to turn off the NaK pumps and the nitrogen supply were placed 
on the main control panel and adjacent to two doors. A future remove valve 
for the cooling-water supply has been provided for In the control panel. 

11.2 NaK-NItrogen Generator Experiments 

The significant parameters for all of the experimental runs for the 
period from August 1975 through September 1976 are summarized In Appendix A. 

The runs were chosen to give operating experience with the revised facility, 
determine the limits of the facility, and allow comparison of measurements with 
the original and revised facilities to see If there were any changes. The data, 
as discussed In Section III, showed a significant Improvement In generator 
performance with the revised facility. 

11.3 Foams 

The slip loss due to the relative phase velocities Is presently one of the 
major loss mechanisms In two-phase liquid-metal MHD generators, especially at 
the high void fractions required for efficient cycle operation. One method of 
reducing this loss Is to produce a more-stable foam flow by means of surface-active 
agents. The potential gains In generator and cycle performances make this 
approach very attractive. Thus, a program to evaluate surface-active agents 
has been Initiated. The major points are: 
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1. Locate suitable agents for the liquid metals of practical Interest. 

2. Determine the techniques required for adequate separation of the 
phases after the generator. 

3. Perform tests In an MHD generator. 

There Is a lack of data In the literature characterizing the preparation and 
properties of a liquid-metal foam of the desired kind, although the available 
literature on the characterization of aqueous and organic forms has given some 
Insight Into the preparation of a liquid-metal foam. In this literature there 
was frequent reference to the use of surfactant additives, and surface tension 
and viscosity appeared to be the relevant physical properties of interest. 
Therefore, scoping studies focused on the effect of surfactants as liquid-metal 
foam Inducers and the potential for fine particle suspensions to act as foam 
promoters have been Initiated. Viscosity was selected as the measured physical 
property for use In correlation of systems of differing bubble behavior. In 
later studies measurements of surface tension will be Included, 

An experimental facility for studying surfactants In room-temperature 
NaK has been built. Initially, the emphasis will be on visual observations 
In Pyrex vessels of the foaming behavior when nitrogen Is bubbled through. 

A glovebox was opened, cleaned thoroughly, and a new balance and the gas system 
for handling ultra-high-purlty nitrogen were Installed, Safety of the NaK 
operations was reviewed. An apparatus for transferring NaK was designed and 
built. 

While waiting for the high-purity NaK, experiments were done with mercury 
to obtain experience with the facility and to develop a "feel” for foam behavior 
In liquid metals, and encouraging results were obtained. The mercury was held 
In a Pyrex vessel. An upturned tube made of copper or of polyethylene having 
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a 0.054 cm orifice was used as the gas bubbler. Visual observations were made 
of the system behavior and a flow rate of 5 cc/mln was adopted as a standard 
condition. Differences In bubble size were noted when the bubbles were formed from 
the polyethylene tube as compared with those formed by the copper bubbler. 

The size differences most certainly relate, at least In oart, to differences 
In the wetting characteristics of the two materials. This observation suggests 
that the wetting behavior between the liquid metal and the mixer may be an 
Important design consideration In a LMMHD generator. 

A system that was found to show some enhancement of bubble lifetime 
was one having one wt % cadmium in mercury. Motion pictures were made of the 
bubble action (45 frames/second). The most oustandlng enhancement of bubble 
lifetime occurred when cesium at a concentration of 3 wppm was used as the 
surfactant. For this mixture (Cs In Hg) up to seven coexisting bubbles were seen, 
with a group of six showing lifetimes up to 10 seconds after the gas flow was 
stopped. Figures II.2 and II.3, excerpted from the movie. Illustrate the difference 
In bubble behavior In the absence and In the presence of cesium. Higher concentrations 
of cesium were also effective, but not to the same degree. These observations 
demonstrate the potential for developing a liquid metal foam system suitable 
for LMMHD application. 

The Initial experiments with NaK wtll focus on the ability of fine 
particle dispersions to promote bubble ltfetlme. Earlier work had Indicated 
that carbon black dispersed In NaK had apparently altered the surface properties. 
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III. GENERATOR PERFORMANCE 


The Initial runs with the second divergent channel In the revised flow 
facility were made for the same conditions as with the old facility to provide a 
check on the operation of the facility. Some of these results are summarized 
In Section 111,1, A detailed comparison for one set of new and old runs for the 
same operating parameters Is given In Section III.2. The effects of reduced 
back pressure and higher flow rates are discussed In Sections III.3 and III.4, 
respectively. A summary of all of the experimental data for the period from 
August 1975 to September 1976 Is given In Appendix A. 

Ill.1 Repeat Runs at the Same Flow Rate 

The liquid velocity, slip ratio, pressure, local load factor, and void 
fraction versus distance along the generator are plotted In Figs. III.1 through 
III.5 for one set of open-circuit runs. No major change with respect to the 
appropriate runs with the old facility. In 1975, was observed. The void fractions 
are somewhat higher than before, and this Is probably connected with the reduced 
slip and higher bubble parameter (drag coefficient), as discussed In the next 
section. There Is, as before, an unusually-hlgh pressure gradient at the generator 
midsection, although this time It seemed less pronounced, and a pronounced drop 
In the load factor, both at the higher void fractions. 

The same quantities are plotted In Figs. III.6 to III.10 for a generator 
with a load resistance of R L * 0.27 mn. Agatn, the variations are similar to 
the results obtained In July 1975. However, the pressure gradient for the 
hlgh-vold-fractlon runs appeared to be more uniform than previously. There Is 
still a pronounced decrease In the local load factor with Increasing void 
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fraction, probably due to the onset of entrainment of liquid droplets by the 
churning motion of the gas In the generator. The fact that the local slip ratio 
Is less than unity at some point Is not significant, this falls within the 
accuracy of the measurement. 

III.2 Comparison of New and Old Runs 

Sets of data with the old loop (7/7/75) and the revised loop Immediately 
after construction (8/18/76 and 8/24/76) are listed In Table III.1. The efficiency 
as a function of quality was virtually unchanged, with the exception of the 
zero-quality runs (see below). However, as shown by the data for a quality of 
about 0,0058, the power level (pressure difference and voltage) was substantially 
higher (about 35% after correction for unequal magnetic fields) for the runs 
Immediately after construction as compared to the runs with the old loop. 

The runs with the highest power level (August 1976) had a lower measured 
average velocity slip ratio K (ratio of gas to liquid velocity). Implying that 
a more-uniform two-phase mixture existed In the generator, the liquid velocity 
was higher (higher voltage, current, and pressure drop), and the exit void fraction 
(ratio of gas volume to total volume) was larger. Also the fluctuation In the 
generator terminal voltage V L was reduced, amounting to only a few percent of 
the average voltage. (Compare this with the previous runs. 4 ) 

The liquid velocities along the generator calculated from the measured 
void fractions and flow rates for the three runs at about 0.0058 quality are 
compared In Fig. III.11. Note that the differences In either the variations of 
or the average liquid velocities are not sufficient to account for the performance 
differences. An Important difference Is In the average load factor F, Table 111,1. 
This Is higher for the runs with the revised facility. Indicating a higher 
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Table III.l 


Comparison of Generator Data with Old and Revised Loops, 
m^ * 6 kg/s. Back Pressure ■ 0.14 MPa Gauge (20 pslg) 



X 

n 9 

B 

r l 

V L 

Ap 

a 

IT 

V 

Date 


[%] 

m 

[mn] 

[V] 

[0.1 MPa] 



7/7/75 

0 

52.9 

1.175 

0.255 

0.506 

2.67 

0 


1.02 

(old) 

0.00580 

49.9 

1.16 

0.254 

0.741 

2.62 

0.49 

1.34 

0.76 


0.00979 

42.1 

1.165 

0.258 

0.819 

2.66 

0.64 

1.217 

0.61 

8/18/76 

0 

39.4 

1.253 

0.26 

0.470 

3.01 

0 


0.89 

(revised) 

0.00216 

48.1 

1.24 

0.27 

0.665 

3.31 

0.27 

1.272 

0.93 


0.00567 

50.3 

1.235 

0.27 

0.867 

3.75 

0.49 

1.047 

0.82 


0.00852 

46.4 

1.238 

0.27 

0.932 

3.90 

0.59 

1.007 

0.74 


0.0110 

43.1 

1.235 

0.27 

1.01 

4.16 

0.68 

0.876 

0.62 

8/24/76 

0 

42.9 

1.223 

0.24 

0.466 

2.98 

0 


0.91 

(revised) 

0.00580 

49.2 

1.227 

0.247 

0,824 

3.69 

0.48 

1.180 

0.83 


0.00825 

44.7 

1.227 

0,248 

0.885 

3.90 

0,55 

1.187 

0.77 


0.01118 

40.6 

1,222 

0.247 

0.937 

4.00 

0.62 

1.18 

0.69 


0.0128 

43.8 

1.225 

0.246 

1.003 

3.99 

0.54 

1.74 

0.89 
















"effective" liquid velocity-electrical conductivity combination for the new runs. 
This result Is expected from the reduced fluctuations. 

It Is suspected that the Improvement In generator performance for the 
tests with the revised loop was due to changes In the surface properties of the 
NaK. Putting the NaK Into the revised loop probably caused a considerable Increase 
In the Impurities In the NaK, primarily Iron oxide and Iron powder from the carbon 
steel. These Impurities could have cuased some surface activity in the NaK, 
and thus Increased the stability of the two-phase mixture and decreased the size 
of the gas bubbles produced by the mixer. 

The zero-quality runs require special comment because the reasons proposed 
above should have no effect. Comparing the July 1975 and August 1976 runs, the 
terminal voltage and efficiency were higher for the former and the pressure drop 
was lower. The measured void fraction profiles Indicate that flow separation 
occurred In at least some of the old zero-quality runs, Including the July 1975 
case. Comparing expected behaviors with and without flow separation, for the 
latter the liquid velocity should be less at the exit end of the generator, the 
circulating currents should be higher, and thus the pressure gradient In the 
upstream part of the generator should be higher. (In the pure liquid case the 
liquid velocity decreases due to the diverging channel, so that the liquid Is 
pimped In the downstream part of the generator.) This predicted behavior difference 
Is consistent with the data. The suspected higher surface activity could result 
In Improved wetting, and thus less or no flow separation. 

The calculated performance from the generator prediction code Is a terminal 
voltage of 0.477 V and an efficiency of 0,425 at 1,2 T, 6.0 kg/s pure NaK, and 
0.25 mn load resistance. The August 1976 runs yielded values close to the 
calculated values. 
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III.3 Reduced Back Pressure 

The liquid velocity, slip ratio, pressure, local load factor, and void 
fraction along the generator are shown In Figs. III.12 through III.16 for a 
back pressure only slightly above atmospheric pressure In the separator tank. 

In the old facility It was not possible to make runs with separator tank 
pressures below about 0.24 MPa absolute (35 psla) due to the high pressure drop 
between the separator tank and the suction side of the pump. With the Improved 
piping system, the lower limit on separator tank pressure Is that It be high 
enough to push the nitrogen through the six-inch line to the atmosphere. 

A reduced generator back pressure means that higher volumetric flow rates 
and velocities are obtained In the generator with the same mass flow rate. 

Also, lowering the back pressure automatically increases the pressure ratio 
across the generator. This In turn means that the liquid velocity will be 
increasing along this generator since It was designed for a lower pressure ratio, 
and this should result In reduced generator performance. 

The local load factor seems to show an Increasing entrainment of liquid 
In the gas along the generator, Fig. III.15. The void fraction and slip ratio 
become quite large at the generator exit. The expected Increase in the liquid 
velocity along the generator Is observed for all of the non-zero quality runs. 

Fig. III.12. 

A comparison of data for normal and reduced generator back pressures in 
Appendix A or Figs. III.6 and 111,12 shows that the terminal voltage V L is 
reduced for the lower-back-pressure runs even though the average liquid velocity 
Is higher. The generator efficiencies for both back pressures are shown In 
Fig. III.17. The efficiency Is lower for the lower back pressure, probably 
because of the varying velocity along the generator and thus Increased circulating 
currents. 
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The ability to go to lower pressure at the generator exit will be useful 
In the future In obtaining larger pressure ratios across the generator, l.e., 
values closer to those anticipated for a commercial generator. Of course, to 
obtain high generator efficiency for higher pressure ratios a generator taper 
designed for the specific pressure ratio Is required. 

III.4 Higher Flow Rates 

The local Darameters for a mass flow rate of 11.9 kg/s are shown In 
Figs. III.IS to III,22. The runs were made at a reduced magnetic flux density, 
0.6 T, so as not to exceed the pressure rating of the channel. The attainable 
quality was limited by some NaK carryover with the exhaust nitrogen, and this 
necessitated Installation of a secondary separator. Section II.1. Note the 
reduced slip ratios and higher load factors that resulted for the higher velocity. 
These are Indicative of a potential Improvement In generator performance with 
higher velocities. 

Generator efficiencies at nominal flow rates of 6 and 12 kg/s with the 
other parameters the same are compared In Fig. III.23. Note that there appears 
to be a tendency for the efficiency to Improve more with Increasing quality at 
the higher velocity than at the lower velocity. (There Is some question about 
the highest quality run at the higher flow rate. Further runs were not made 
before the end of the period covered by this report at this quality due to NaK 
carryover.) 
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Figure III.20 Pressure Along Generator, 11.91 kq/s 
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DISTANCE ALONG GENERATOR (m) 

Figure III.21 Local Load Factor Alonq Generator, * 11.91 kg/s 
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IV. ANALYTICAL STUDIES OF ANNULAR GENERATORS 


IV.1 Introduction 

Rectangular MHD conduction machine geometries for generators and liquid- 
metal pumps have been studied almost exclusively. Helical geometries are 
suitable for low-volume-flow-rate high-pressure difference pumps, but not for 
MHD generators. Annular geometries are common for Induction pumps and generators, 
but there are no known references on annular conduction machines. Since the 
annular geometry has several advantages, the following analytical study was 
undertaken to evaluate Its practicality. 

The annular generator geometry considered consists of two concentric, 
homogeneous cylindrical walls of different radii with the conducting fluid flowing 
In the axial direction within the annulus. Fig. IV.1. Advantages of this geometry 
are: 

1. The cylindrical walls are the electrodes, eliminating the need for 
Insulating walls and the resulting fabrication problems. Without the Insulating 
walls there are no accompanying viscous boundary layer and current shunting 
effects as found In rectangular geometries. 

2. Homogeneous cylindrical shells are Ideal pressure vessels. 

3. The circumferentially-directed magnetic field Is entirely contained 
within the generator, stray fields are not a problem, and the field Is produced 
by a superconducting or ordinary torroldal winding without any magnetic Iron. 
Dlsdavantages of this geometry are: 

1. The electrical terminal properties are Inherently low voltage and high 
current because of the relatlvely-short electrode spacing and large cross-sectional 
area for current flow. 
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EXCITING TOROIDAL WINDING, 
NI e , c AMPERE-TURNS 







2. The electromagnetic force density has an Inverse-square radial dependence 


which may result In a very nonuniform velocity profile. 

3. The magnetic field generated by the Induced currents (armature reaction) 
may cause large variations In the power density along the machine, and there Is 

no apparent method to compensate for this without destroying the first advantage. 

4. The toroidal magnet winding has to be closed through the end regions 
of the annular duct, so that the wires must cross the flow area. 

The first three disadvantages are analyzed In the following material by 
means of three simplified models. After establishing the relevant equations 
and assumDtlons, a model with a constant velocity and no armature reaction Is 
considered to establish practical regions of operation In Section IV.3. A 
laminar velocity profile without armature reaction. Section IV.4, Is used to 
study the effect of the Inverse-square radial dependence of the electromagnetic 
force density. The flow In a practical generator could be either laminar or 
turbulent, depending on the relative lamlnarlzlng effect of the magnetic field, 
but the laminar solution yields Insight Into what might happen even for turbulent 
flow. Finally, the axial dependences of the magnetic flux and current densities 
for a constant fluid velocity are considered In Section IV.5. Numerical results 
obtained from an annular-geometry computer code are Included for all three 
models. The primary objective of this study Is generator operation, but the 
solutions apply equally to oumo operation and numerical oump cases are also 
discussed. Only pure-liquid machines are considered, but the conclusions are 
valid for two-phase liquid-metal and plasma machines as well. The velocity- 
profile solution also has Important Implications for annular Induction pumps 
and generators. 
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IV.2 Equations and Assumptions 

All models considered are qoverned by the same general equations of 
fluid mechanics and electromagnetics. The assumptions made to solve the 
equations are: 

1. Steady-state conditions, 

2. Incompressible, Newtonian, non-polarlzable fluid, l.e., pure 
liquid, 

3. No fluid sources or sinks In the generator, 

4. Negligible gravitational effects, 

5. One-dimensional axial flow with circular symmetry about the axis 
of the cylindrical walls, and 

6. The magnetic permeability Is that of free space everywhere. 

With these assumptions, the basic relevant equations are 


V X § * v 0 J, 

(IV.1) 

3 a o(t + U X ft). 

(IV.2) 

v x f = 0, 

(IV.3) 

V • $ - 0, 

(IV.4) 

v • u - 0, 

(IV.5) 

p(u • v)u = -vp + pv 2 u +1x5; 

(IV.6) 


where the variables are defined In the nomenclature. 

In cylindrical coordinates with B 0 , u z , J r , and E r as the only vector 
components the equations become 

IF (re 9 ) . 0, (IV.7) 



(IV.8) 


(IV. 
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where NI gxc Is the ampere-turns of the exciting torroldal winding. Fig. IV.1. 

The electrical terminal properties, the electrode or load current 1^ 
and the electrode voltage difference or load voltage V L , are given by 



-2irr | J r dz 
0 

o 


(IV.16) 

(IV.17) 
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where L Is the length of the machine, and r. and r are the radii of the Inner 

I o 

and outer cylinders (electrodes). 

The general equations for the powers are collected here for reference. 
The ohmic loss in the fluid Is 


P 


r 


L 



0 


rdrdedz. 


(IV.18) 


The electrical output power with end loss neglected (negative for a pump) Is 
p e = Vl* (IV.19) 

The mechanical Input power (negative for a oump) is 

P m * Q(Ap em + Ap v >’ (IV.20) 

where Q Is the volumetric flow rate of the working fluid, Ap^ Is the electro¬ 
magnetic pressure difference, 

fpLu 2 

Ap v a "2I5 r~ ( lv - 21 ) 

h 

Is the (non-MHD) viscous pressure difference In terms of the Darcy friction factor 
f, and ■ 2(r Q - r^) Is the hydraulic diameter. The efficiencies for generator 
and pump operation are defined as 


n g 5 


_ P m 

» = 'PT 


(IV.22) 


It Is convenient to use the normalized or non-dimensional radius 
R = r/r 1 (IV.23) 

In the analysis. At the outer electrode, 
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(IV.23a) 





r 


IV.3 Constant Velocity with no Armature Reaction 

Constant fluid velocity with no armature reaction (the Ideal case) 

Is assumed to establish the parameters required for practical applications. The 
required equations are determined, and numerical examples calculated for both 
generator and pump operation. 

IV.3.1 Analysis 

The load or terminal current Is 

I L » -2nrLJ r (IV.24) 

from Eq. IV.16. The terminal voltage Is calculated from Eq. IV.17, where 
E r Is obtained by substituting 1^ and B 0 from Eq. IV.15 Into Eq. IV.9 and 
solving for E r . The result Is 


' u z “."m 

K 


2ttLo 


lnR„ 

o 


(IV.25) 


In terms of the normalized radius. This can be Interpreted as the equivalent 
circuit of Fig. IV.2, which has the equation 


V L * V oc - \*V 
The open-circuit voltage Is 


(IV.26) 


(IV.27) 


and the Internal resistance Is 

R 1 a lnR 0 /(27rLo). (IV.28) 

In terms of the magnetic flux density at the Inner radius, Bj, V QC can be 
written as 

V * ".WV (IV - 27,) 
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(IV.29) 


The ohmic power loss In the fluid, from Eq. IV.18, Is 

P r " I L 2 V 

The electrical output power, Eq. IV.19, Is written In terms of V L as 
P e * oLV L (y u 2 NI exc - 2n V L /lnR Q ) (IV.30) 

using Eq. IV.25 or Eq. IV.26 to eliminate I L . In terms of the load factor, F, 
defined as 


F -w 

Eqs. IV.29 and IV.30 can be written as 


P r * (Voc^Kl - F ) 2 

and 


(IV.31) 

(IV.29a) 


P e “ ( V oc 2/R l) F( 1 “ F) * (IV.30a) 

The calculation of P m Is more complicated because the electromagnetic pressure 
gradient Is not constant over the channel cross section. To calculate Ap^, 
the Lorentz force density 3 x § Is Integrated over the total channel volume, 
yielding 


F «. ’ -” N W W 2 "' < IV - 32 > 

The minus sign means that F^ acts In the -z direction or opposes the fluid motion. 
Then, 

-F u NI I, InR 

3 —r 1 ■ ——-- (IV.33) 


2*2 r^ (R q 2 - 1) 


Is the average difference. The mechanical Input power Is calculated from 
Eqs. IV.20, IV.21, and IV.32, and the efficiency from Eq. IV.22. If the 
viscous oressure difference Is negligible, then 

P m " <V 0C 2/ Ri>(l - F) (IV.34) 

and 

n g » F (IV.35) 

as In rectangular MHO machines. 
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IV.3.2 Numerical Examples 

The most significant parameter Is the obtainable terminal voltage, as 
must be large enough for either efficient dc to ac Inversion or homopolar 
motors. Since the load factor F must be close to unity for efficient generator 
or pump operation, Eq. IV.35, V oc Is given In Table IV.1 for values of r^, R Q , 
and u z . Note that R Q Is limited to be much less than 1.1 at reasonable power 
densities to prevent the possibility of flow reversal, as shown In Section IV.4. 
The velocity Is limited to around 25 to 50 m/s by viscous loss. None of the 
cases shown yields adequate voltage, even at = 6 T and r^ = 5 m. Higher 
voltages are obtainable with several generators connected electrically In series, 
except that the output power per generator (see Table IV.2) may be too large 
to use several generators In one power plant. 

The powers are shown In Tables IV.2 as a function of the load factor F 
for one of the cases from Table IV.1. For an acceptable efficiency with all 
losses Included (end losses are Ignored here) F will most likely be between 
0.9 and 0.95, leading to a generator power on the order of 500 MW per meter 
of length and V L around 7 V at r^ 5 1 m and = 3 T. Since the power scales 
as V QC 2 (or r^ 2 or B^ 2 ) at constant R Q and F, Eq. IV.30a, Increasing V QC 
results In Increased P g and no possibility of a series connection to raise the 
useable voltage. Flndlnq a good compromise between generator size, power level, 
and terminal voltage Is very difficult for a rectangular generator, and still 
more dlfflcult.lf not Impossible, for larqe annular generators. 

There are no applications for electromagnetic pumps of the power levels 
shown In Table IV.2. Thus, acceptable voltage levels may be unattainable for 
annular conduction pumps since the size and voltage must be reduced from the 
values of Table IV.2. The low voltages shown \i Table IV.2 are not surprising, 
they result from the relatlvely-short electrode spacing mentioned In Section IV.1 
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Table IV.2 




Powers per Meter Length 

Fluid 3 NaK 
B< 3 3 T 
25 m/s 
1.5 x 


7.15 V 


for r^ 
10 7 AT 


1 m, R Q = 1.1, and Constant Velocity 

p 3 867 kg/m 3 
<j « 2.62 x 10 6 mhos/m 
D h 3 0,2 
f"« 0.022 
= 82,3 



F 

V L 

[V] 

! l 

[MA] 

P e 

[MW] 

P m 

[MW] 

Useful Power 
Density 
[MW/m 3 ] 

n 9 

m 

Gen. 

0.8 

5.72 

247 

1410 

1770 

2140 

79.9 

Gen. 

0.95 

6.79 

61.8 

420 

442 

636 

94.9 

Pump 

1.1 

7.86 

-123 

-971 

-882 

1340 

90.9 

Pump 

1.3 

9.29 

-370 

-3440 

-2650 

4010 

76.9 
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In this section the effect of the radial dependence of the electromagnetic 


force density Is analyzed. Fully-developed laminar flow Is treated because 
analytical solutions can be obtained. The laminar-flow solution provides a 
guide for the behavior of turbulent flows. 

IV.4.1 Analysis 

First the velocity profile u z (r) must be determined. The electromagnetic 
force density O r B 0 can be written as 


% NI exc ! L . -1 B 1 *L 
(2nr) 2 L 2*r 2 L 


(IV.36) 


from Eqs. IV.15 and IV.24. The pressure gradient can be expressed as 3p/3z = -Ap/L, 
where Ap equals the Inlet minus the exit pressures. Note that Ap Is Independent 
of r and z for fully-developed flow and Is positive for generator operation. 
Substituting Into Eq. IV.13 and solving for u 2 (r) with the boundary conditions 
that u z = 0 at the walls yields 


u (R) * [-In R ln(R /R)] - 


R 2 - 1 


In R_ 


- In R - (R 2 - 1) 


(IV.37) 


In terms of the normalized radius R. Note that this Is not a Hartmann profile, 
and there Is no parameter like the Hartmann number governing the shape of the 
velocity profile. 

Next V L Is calculated using the same procedure as In Section IV.3.1 
except that u z (R) Is given by Eq. IV.37. The result. In terms of the equivalent 
circuit of Fig. IV.2 or Eq. IV.26, Is 






(IV.33) 

(IV.39) 


Note that R^ Is not the physical internal resistance In R q /2itLct. The second 
term, which is independent of o, comes from the specified constant pressure 
difference. With Ap fixed, the average fluid velocity or the volume flow rate 
Q decreases as 1^ Increases because of the increased (adverse) electromagnetic 
pressure gradient. In Section IV.3.1, on the other hand, u z (or Q) was specified 
and Ap had to increase with increasing 1^. 

The ohmic power loss Is still qiven by Eq. IV.29 since the second term 
of Eq. IV.39 does not represent a physical resistance. The electrical output 
power P g = V L I^ is expressed as a function of Ap and I L using Eqs. IV.26, IV.38, 
and IV.39. The mechanical input power is calculated from P m = Qap, where the 
volumetric flow rate is 


Q 



Integrating Eq. IV.40 using Eq. IV.37 results in 


(IV.40) 



(IV.41) 


The viscous power dissipation can be calculated from 
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r o 3U 

P tf * 2 itL v ( -r^- )2 rdr (IV.42) 

r, 

using Eq. IV.37, or from the other powers and conservation of energy. 

IV.4.2 Numerical Examples 

The velocity profiles for r^ s 3 m, R Q * 1.01, = 3 T, V L = 4 V, NaK, 

and three pressure gradients are shown In Fig. IV.3. The powers for these 
parameters and six pressure gradients. Including the three plotted In Fig. IV.3, 
are given In Table IV.3. Increasing pressure gradient corresponds to decreasing 
load factor F, as shown In Table IV.3. Note that for small pressure gradients 
(F close to unity) the profile resembles the non-electromagnetlc annular 
profile, the viscous and ohmic losses are small, and the power density Is low. 

At a power density of about 30 MW/m 3 flow reversal occurs (Fig. IV.3) and the 
viscous power loss Is substantial. Although the efficiency Is very high, the 
load factor Is so close to unity that the power density is relatively low for 
practical MHD generators and the end losses (not Included) would be serious. 

A similar behavior Is observed for R Q ■ 1.1, Fig. IV.4, except that the 
flow reversal occurs at a lower pressure gradient or power density because of 
the larger variation In the electromagnetic pressure gradient across the channel 
cross section. 

Note that the behavior of the annular generator Is different from the 
convertlonal rectangular generator In that flow reversal occurs and the viscous 
losses Increase sharply with generator loading (decreasing F). Flow reversal 
occurs for all cases except very low power densities so that the viscous 
contribution to the total pressure qradlent can balance out the radial variation 
of the electromagnetic pressure gradient. In the rectangular generator 
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Figure IV.3 Velocity Versus Radius as a Function of the Pressure Gradient, 









VELOCITY 



Flqure IV.4 Velocity Versus Radius as a Function of Pressure Gradient, R 










Table IV.3 


Powers per Meter Length and Efficiencies as a Function of the 
Pressure Gradient for r 4 - 3.0 m, R a 1.01, 

Bj a 3 T, V|_ - 4 V, and NaK 


AP/L 

[N/m 3 ] 

F 

*9 

m 

Useful Power 
Density 

[MW/m3] 

P e 

[MW] 

P n, 

[MW] 

P v 

[MW] 

p r 

[MW] 

6.9 x 10 2 

1.0000 

27.0 

0.00832 

0.00473 

0.0175 

0.0128 

8.44 x 10' 10 

6.9 x 10 3 

1.0000 

92.7 

0.286 

0.162 

0.175 

0.0128 

9.96 x 10" 7 

6.9 x 10 4 

0.9999 

99.2 

3.06 

1.74 

1.75 

0.0144 

1.14 x 10' 4 

6.9 x 10 5 

0.9993 

99.0 

30.8 

17.5 

17.7 

0.171 

1.16 x 10‘ 2 

6.9 x 10 6 

0.9934 

91.2 

308. 

175. 

192. 

15.8 

1.16 

6.9 x 10 7 

0.9379 

50.8 

3080. 

1750. 

3450. 

1580. 

116. 
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viscous power loss Is almost constant, except for the slight Increase In average 
velocity needed to hold the load voltage constant as F Is decreased. In the 
cases considered for the annular generator the viscous loss Is always much larger 
than the ohmic loss, and the efficiency decrease at higher power density 
(lower F) Is caused by the sharp Increase In the viscous loss with the flow 
reversal and the resulting large shear forces. 

The velocity profiles for two sets of pump cases, R Q = 1.01 and 1.1, are 
shown In Figs. IV.5 and IV.6, respectively. The corresponding powers and 
efficiencies are given In Tables IV.4 and IV.5. In both cases flow reversal 
occurs, although at a much higher pressure gradient or power density for the 
R q ■ 1.01 case because of the smaller change In the electromagnetic force density 
across the channel. Note that for even higher pressure gradients or power 
densities the devices become dlsslpators, absorbing both electrical and mechanical 
Input powers and dissipating them as ohmic and viscous losses. (A similar 
effect will occur with any pump. A pump has a given pressure difference with 
zero flow; If a higher pressure difference Is Imposed, then reversal of the net 
flow will occur.) 

The voltage levels for the pump cases are very low because of the short 
electrode spacing. Both cases considered have the same average electric-field 
strength. The low voltages mean that the device Is probably not practical. 

The results obtained here are only valid for laminar flows. However, 
similar effects are expected with turbulent flows. Including flow reversal 
at high power densities. 
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VELOCITY U z (m/s) 



Flqure IV.5. Velocity versus Radius as a Function of the Pressure Gradient 
for a Pump, R ■ 1.01 













VELOCITY U z (m/s) 



Flqure IV.6 Velocity versus Radius as a Function of the Pressure Gradient 
for a PumD. R. * 1.1 











Table IV.4 

Powers per Meter Length and Efficiencies as a Function of the 
Pressure Gradient for a Pump, r* ■ 0.5 m, R 0 » 1,01, 

B 1 ■ 2 T, V L - 0.05 V, and NaK 


Ap/L 

[N/n 3 ] 

F 

n P 

[%] 

Useful Power 
Density 
[kW/m 3 J 

-Pe 

[KH] 

"Pm 

CMC 

p y 

Uffl 

P r 

m 

-6.9 x 10 2 

1.000 

25 

3.466 

0.2166 

0,05473 

0.1619 

1.13 x 10" 5 

-6.9 x 10 3 

1.000 

77 

34.66 

0,7092 

0.5472 

0.1618 

1.22 x 10" 4 

-6.9 x 10 4 

1.001 

97 

346.1 

5.635 

5.465 

0.1626 

7.6 x 10' 3 

-6.9 x 10 5 

1.013 

98 

3,413 

54.89 

53.89 

0.2798 

0.7286 

-6.9 x 10 6 

1.153 

85 

29,300 

547.5 

462.7 

12.34 

72.47 

-6.9 x 10 7 

- 

- 

- 

5,473 

-2,992 

1,222 

7,243 


Table IV.5 

Powers per Meter Length and Efficiencies as a Function of the 
Pressure Gradient for a Pump, r+ « 0.5, R 0 ■ 1.1, 

B| * 2 T, V L a 0.5 V, and NaK 


AP/L 

[N/m 3 ] 

F 

o n Useful Power 
p Density 

L%1 [kW/m3j 

" P e 

[kW] 

" P m 

[kW] 

P v 

CkVO 

P r 

[kW] 

-6.9 x 10 2 

1.000 

95 

3.73 

0.6151 

0.5849 

0.0301 

8.7 x 10‘ 

-6.9 x 10 3 

1.000 

80 

36.30 

5.985 

4.796 

1.189 

8.3 x 10‘ 

-6.9 x 10 4 

1.001 

- 

- 

59.68 

-57.44 

117.0 

0.0825 
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IV.5 Armature Reaction with Constant Velocity 

The effect of the axial variation of the magnetic flux density and the 
(electric) current density caused by the current In the fluid Is considered 
here. The three different connections to complete the load circuit, shown In 
Fig. IV.7, are: 

1. At the downstream end, 

2. At the uDStream end, and 

3. At both ends. 

The three connections, their different effects on the magnetic field, and the 
resulting performances are discussed In the following material. 

IV.5.1 Downstream Connection 

If the load Is connected to the downstream ends of the electrodes, then 
the z-dlrected load current flowing along the Inner electrode Is In the same 
direction as the exciting current for a generator. Fig. IV.7. Thus, the 
Induced magnetic field due to the load current reinforces the applied field, 
and the total field Is every place greater than the aonlled field (except at the 
generator entrance where the Induced field Is zero). 

The apolled field Is riven by Eq. IV.15, and the total field Is 

u n N W r 

B e (r,z)- Vif? + u oJ J r (r ’ z ' )dz ' (IV.43) 

Q 

from Eq. IV.14 and conservation of charge. The electric field E f must vary 
as 1/r from conservation of charge or Eq. IV.8,and Eq. IV.9; and In terms of V L , 








(c) BOTH ENDS CONNECTION 


Figure IV.7 Load Connections for Armature Reaction Calculations 










from Eq. IV.17. Substituting Into Ohms law, Eq. IV.9, gives 

% r mi. 


-a V. * 

•\>- 2 > ■ fdt ♦ r 


•1 


■Ur.z')dz> . 


(IV.45) 


R m = V u z L (IV.46) 

Is the magnetic Reynolds number. Solving Eq. IV.45 yields 


J„(r,z) . 


m exc p L 
TirrE rTrTT 


\ Z ' L . 


(IV.47) 


and, from Eqs. IV.43 and IV.16, 


NI 


B 9 < r - 2 > * 2/L + a ~ rlnV N ' <"- 48 > 


“a 2 / 1 -. 


t“ NI exc ’ ~uTTnH* (e " 1) ‘ 


(IV.49) 


From comparing Eq. IV.49 with the equivalent circuit, Fig. IV.2, V QC Is given 
by Eq. IV.27 for the Ideal case, constant velocity with no armature reaction, 
since there Is no armature reaction for open circuit. The Internal resistance 
for this case. 


R_ In R 

R » m _ 0 

R 1 ft 

2*Lo(e ™ - 1) 


(IV.50) 


Is always less than the value for the Ideal case, Eq. IV.28. The reduced value 
for Rj occurs because the total magnetic flux In the generator Increases with 
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loading, and this acts to compensate for the Increased ohmic voltage 
drop. 

The Lorentz force Integrated over the total channel volume Is 


-v In MNI av J 2 (l - F) R m 

-2- fa-*™ -[2F(e - 1) + (1 - F)(e - 1)]. 

(IV.51) 


From this plus the viscous pressure drop, Eq. IV.21, the mechanical Input power 
can be calculated. The ohmic power loss In the fluid Is 


% * ( V - V 2 (e2Rm -')- < IV - 52 > 

The remaining powers and efficiencies can be calculated from the equations 
In Section IV.2. 

IV.5.2 Upstream Connection 

If the load Is connected to the upstream ends of the electrodes, then the 
z-dlrected load current flowing along the Inner electrode Is In the opposite 
direction to the exciting current, Flq. IV.7. For this case the Induced magnetic 
field opposes the applied field, the total field Is everyplace less than the 
applied field (except at the generator exit where the Induced field Is zero), 
and the reduction In the total magnetic flux In the generator Is reflected In a 
larger Internal resistance. 

The procedure used to determine the equations for the upstream connection 
Is the same as for the downstream connection. The results are: 
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B.(r.z) 


(IV.43a) 


2 nr 


J (r,z')dz' 


J r (r,z) 
Bfl( r »z) ■ 


‘m exc L 
ZnrL rTFlT 


V 2 ' L -’> 


(IV.47a) 


M, exc A ,2 ' l -’> »l n > (2/L ‘ 

* u ~ ; r l n H. fl ’ e >" 


2irr 


e xc - P 0 u 2 in R 0 
R m R o 


(1 - e"%. 


2irLa(l 


7S’ 


(IV.43a) 
(IV.49a) 

(IV.50a) 


and 

F an ■ fe ln R o (NI exc )2( ' " F >t<' ' F ><' ■ e^S + 2FO - e'S]. 

(IV.51a) 

Note that V QC Is unchanged, and Is always greater than the value for the 
Ideal case, Eq. IV.28. 


IV.5.3 Both Ends Connection 

With the load connected to both ends of the generator. Fig. IV.7, the 
generator Is split Into two regions by the location z * L' where the component 
of the load current flowing along the axis Is zero. The region for z < L' 
behaves like a generator with the upstream connection, the region with z > L' 
behaves like a generator with the downstream connection, and the total properties 
can be determined from the sum of the two regions using the equations from 
Sections IV.5.1 and IV.5.2. For example, the load current Is given by 
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u. In R 


-RJ-'/L 
e m ), 

(IV.49b) 


where Is based on the total length L. 

The location of L' would normally be determined by the voltage drop 
along the Inner electrode. Since the electrode Is assumed to be an equlpotentlal 
here. It Is not possible to calculate a value for L'. 


IV.5.4 Numerical Example 

The powers and efficiencies for one set of parameters are shown In 
Table IV.6 for the three connections, and the magnetic flux densities and 
current densities are shown In Figs. IV.8 and IV.9. For comparison, the values 
with the armature reaction Ignored are also shown. The value of R^ = 4.939 
Is not large, yet the effects are dramatic. In fact, the downstream connection 
results are not realistic, the pressure gradient and magnetic flux density 
both become unrealistically large. (No Iron Is used, so that magnetic saturation 
Is not a factor.) The upstream connection Is realizable, but the power 
density Is low. The choices u z = 15 m/s and o * 2.62 x 10 5 mhos/m (0.1 times 
a of NaK, characteristic of a hlgh-vold-fractlon flow) were made to restrict 
R m to small values. The excitation was chosen to .yield = 3 T with no armature 
reaction. 

Similar effects due to armature reaction occur for rectangular geometries. 
However, It Is easy to comoensate rectangular generators by returning the load 
current through the magnet air gap. For an annular geometry, compensation 
would mean running the load current through the fluid gap between the electrodes, 
and this would eliminate one prime advantage of the annular geometry—the lack 
of Insulating walls between the two electrodes. 
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Powers and Efficiencies With and Without Armature Reaction, 
ri-l.On, Rfl-l.I.L* 1 m, u 2 * 15 m/s, Vi 3 3.903, 

F 3 0.91, o * 2.62 x 10 5 mhos/m, and » 4.939 



Armature 


Connection 


Reaction 

Ignored 

Downstream 

Upstream 

Both Ends, 

I L [MA] 

4.445 

124.7 

0.8936 

2.399 

P e [MW] 

17.92 

502.8 

3.602 

9.673 

P m [MW] 

19.17 

2,759 

3.824 

10.61 

P r [MW] 

1.14 

2,256 

0.1158 

0.835 

P d *[MW/m 3 ] 

27.16 

762.1 

5.46 

14.66 

Ap/L [N/m 3 ] 

1.937 x 10 6 

2.789 x 10 8 

3.865 x 10 5 

1.073 x 10* 

n g [*] 

93.5 

18.2 

94.2 

91.1 

B 1 (z-0) [T] 

3.0 

3.0 

2.821 

2.823 

B< (z-L) [T] 

3.0 

27.94 

3.00 

3.303 


*Pj 3 useful power density 









Bj (T) 



with Arweture Reaction 





(r = Ti) (A/m ) 



Figure IV.9 Current Density at Inner Radius with Armature Reaction 
Versus Axial Location 







V. CONCLUSIONS AND FUTURE WORK 


The ambient-temperature NaK-nltrogen facility has been successfully revised 
to handle higher velocities and pressures, and satisfactory operation has been 
obtained. The Improved control features make the system easier to use. New 
experiments with the second dlverglng-'hannel generator showed that the power 
density was higher, and that there appeared to be ? tendency for the efficiency 
to Improve more with Increasing quality at higher velocities. 

The program to produce and evaluate liquid-metal foams was Initiated. 

Equipment was developed, screening tests made with mercury while waiting for 
the high-purity NaK, and a olan established. The future program Is described 
briefly below. 

The evaluation of an annular generator geometry was completed. Although 
this geometry appears advantageous because there Is no need for Insulating walls, 
the cylindrical (electrode) shells are Ideal pressure vessels, and the magnetic 
field Is contained entirely within the generator; there are also serious limitations. 
The limitations are; 

1. The geometry Is Inherently suited to low voltage high-current operation, 
not to the higher terminal voltages needed as Input to acyclic motors or to 
efficient and cheap dc to ac converters. 

2. Flow reversal and Increased viscous loss can occur at the power densities 
desired for power systems. 

3. Armature reaction effects are significant and compensation methods 
appear Impractical. 

Although these limitations were known before the analysis described In Section IV 
was performed, the analysis was necessary to quantify the effects. Based on this 
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analysis, this annular geometry does not appear attractive for either generator or 
pump operation. 

The desired next steps In the liquid-metal MHD program are described 
below. 

V.l Higher Velocities and Pressures 

The motivations for tests at higher liquid velocities are to reduce the 
slip loss due to the relative ohase velocities, as experimental research has 
demonstrated that this loss can be decreased by Increasing the liquid velocity, 3 
and to test the generator at velocities typical of practical generators. In 
the experiments to date except for a few recent runs, the liquid velocity has 
been limited to about 12 m/s. The revised facility allows generator operation 
at higher liquid velocities or flow rates and higher pressures. Testing of the 
existing diverging channel will continue at liquid flow rates up to about 
10.9 kg/s (200 gpm), versus the previous limit of about 6 kg/s (110 gpm). This 
generator was designed for lower velocities and has an operating pressure range 
up to about 0.72 MPa absolute (90 pslg), so that operation at higher velocities 
will require a reduced magnetic field strength. 

A new generator will be designed and built to utilize the full loop 
capability of 10.9 kg/s and 1.48 MPa absolute (200 pslg). This will allow 
high velocities with high field strengths, and variation of the generator 
pressure range with a maximum generator exit pressure of 1.14 MPa absolute 
(150 pslg). This generator will have a different taper to fit the different 
design pressure ratio. Also, the maximum void fraction In the generator should 
be Increased—this will be aided by a reduced velocity slip ratio. 
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V.2 Foams 

The Initial steps of the program to evaluate foams In MHD generators were 
described In Section II.3. In the next year attention will be focussed on 
Identifying surface-active agents suitable for NaK, and plans should be developed 
for testing NaK foams In a generator. Bubbling tests with NaK and various 
additives will be part of the screening program, along with measurements of 
surface tension, viscosity, and wetting oropertles. NaK was chosen for the 
Initial tests because of extensive experience with NaK In generators, and thus 
the ability to see the change In generator performance due to the use of a 
foam In minimum time and for minimal cost. If the NaK test results are positive, 
then effort will be directed toward finding foams suitable for use In practical 
MHD systems. Methods of separating the gas and liquid phases after the MHD 
generator for a foam must also be considered. 

V.3 Improved Local Diagnostics 

Improved local diagnostics are needed to obtain a better understanding of 
two-phase flows In a magnetic field. The data to be measured include direct 
local values of liquid and gas velocities, void fractions, bubble diameters, 
pressures, and temperatures. Measurement techniques to be evaluated Include: 

(1) conductivity orobes, both single and multiple, for void fractions, bubble 
sizes, and bubble velocities; (2) hot-film probes for liquid and gas velocities, 
void fractions, and bubble sizes; and (3) microthermocouples for gas and liquid 
temperatures. These techniques have been used successfully for air-water two-phase 
flows, but have not been tried with liquid metals or magnetic fields. Special 
electronics and data-processlng techniques may be required. 
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This local data Is needed to develop and evaluate Improved analytical models 
for two-phase flows. Section V.4. Thus, a program to develop the measurement 
techniques should be Initiated. First tests should use the air-water facility 
described In Reference 4 before Incorporation Into the ambient-temperature 
NaK-nltrogen loop. 

V.4 Analytical Models for the Two-Phase Flow In the Generator 

Further analytical studies on the two-Dhase flow In the MHD generator should 
be carried out In parallel with the above efforts. The experimental data, 
especially that from the Improved local diagnostics as It becomes available, 
w111 be used In developing and evaluating these models. The primary goals are 
a better understanding of and model for the two-phase flow at high void fractions 
In a practical generator. The conditions for transitions between bubbly, 
churn-turbulent, slug, and annular flows at high void fractions may be explored. 

Of particular Importance are the effects on generator performance of velocity, 
void fraction, and conductivity distributions across the channel, l.e., parallel 
to either the applied magnetic field or the current. These have not been 
Included In any model to date except for a simple boundary-layer analysis. 3 

V.5 Imperfect Compensation 

Most analyses of two-phase MHD generators have assumed perfect compensation, 
l.e., the generator current Is returned through the air gap In the magnetic 
circuit In such a fashion that the applied magnetic field Is not changed by the 
generator current. This may not be true In practice because of end currents, 
non-uniform current distributions In the generator and/or compensating conductors. 
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and compensating conductors that are longer than the electrodes as In the 
present generator. The development of an analytical model for the generator 
with Imperfect compensation has been Initiated. Future plans Include experi¬ 
mental verification of the model. 
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APPENDIX A 

SUMMARY OF EXPERIMENTAL GENERATOR DATA 
AUGUST 1975 TO SEPTEMBER 1976 


Date 

Run 

Number 

Mass Flow Rate 
NaK R 2 

[k"js] [kg?s] 

Magnetic 

Flux 

Density 

B 

[T] 

Load 

Resistance 

r L 

[mn] 

Load 

Voltage 

c?i 

Pressure 
Inlet Exit 

p 1n R out 

[psla] [psla] 

Power 

Output 

P e 

[kW] 

Average 

Void 

Fraction 

sr 

8-12-76 

1 

6.025 

0 

1.245 

- 

.509 

63.3 

37.5 

0 

0 


2 

6.115 

.0125 

1.246 

- 

.717 

57.7 

34.2 

0 

.277 


3 

6.025 

.03175 

1.244 

- 

.932 

57.1 

35.2 

0 

.480 


4 

6.025 

.0507 

1.245 

- 

1.102 

59.5 

36.8 

0 

.558 


5 

6.070 

.0655 

1.245 

- 

1.257 

60.7 

36.0 

0 

.652 

8-16-76 

1 

12.140 

0 

1.245 

- 

.996 

88.7 

31.0 

0 

0 


2 

12.095 

.0135 

1.235 

- 

1.316 

74.3 

26.4 

0 

.184 


3 

11.560 

.0359 

1.234 

- 

1.784 

78.4 

31.0 

0 

.450 

8-18-76 

1 

6.160 

.0351 

1.235 

.2706 

.867 

88.4 

34.0 

2.78 

.495 


2 

6.003 

.0516 

1.238 

.271 

.932 

92.6 

36.3 

3.21 

.587 


3 

6.025 

.0672 

1.235 

.271 

1.010 

93.5 

33.7 

3.76 

.678 


4 

6.025 

.0130 

1.240 

.2703 

.665 

80.3 

31.5 

1.63 

.273 


5 

6.025 

0 

1.253 

.264 

.470 

81.1 

35.0 

.837 

0 


6 

6.025 

.0654 

1.235 

.270 

.931 

77.0 

18.5 

3.21 

.672 

8-19-76 

1 

6.038 

0 

1.248 

.2644 

.476 

65.0 

18.3 

.856 

0 


2 

6.025 

.01396 

1.233 

.2709 

.852 

70.7 

15.5 

2.70 

.524 


3 

6.025 

.0356 

1.230 

.2704 

.725 

65.4 

16.4 

1.94 

.388 


4 

6.025 

0 

1.25 

- 

.492 

45.1 

18.1 

0 

0 


5 

11.825 

0 

1.265 

- 

.982 

64.3 

15.5 

0 

0 


6 

6.070 

.01231 

1.236 

- 

.827 

39.2 

17.5 

0 

.435 


7 

6.025 

.03078 

1.235 


1.102 

36.8 

21.7 

0 

.608 
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II 


Run 

Date Number 

Mass Flow Rate 
NaK *2 

Magnetic 

Flux 

Density 

B 

[T] 

Load 

Resistance 

r L 

[mh] 

Load 

Voltage 

eft 

Pressure 
Inlet Exit 

p 1n p out 

[psla] [psla] 

Power 

Output 

P e 

[MO 

Average 

Void 

Fraction 

ar 

9-17-76 1 

6.151 

0 

.6 

.250 

.232 

47.0 

36.1 

.215 

0 

2 

6.002 

.01261 

.6 

.255 

.344 

49.5 

36.4 

.464 

.317 

3 

5.98 

.03318 

.595 

.258 

.451 

52.0 

36.8 

.788 

.520 

4 

6.07 

.06330 

.595 

.259 

.539 

53.8 

35.1 

1.12 

.610 

5 

11.915 

.07023 

.588 

.373 

1.068 

64.2 

32.0 

3.06 

.541 

9-20-76 1 

5.980 

0 

.599 

- 

.238 

43.9 

36.6 

0. 

0 

2 

6.025 

.0135 

.613 

- 

.403 

41.8 

35.4 

0 

.359 

3 

6.025 

.0348 

.598 

» 

.543 

41.7 

35.8 

0 

.585 

4 

6.025 

.0502 

.598 

« 

.638 

42.9 

34.4 

0 

.637 

5 

5.98 

.06378 

.602 

. 

.695 

43.7 

35.5 

0 

.726 
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